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ABSTRACT 

In a 6h Ha exposure of the N-W region of the cluster of galaxies A1367 we discovered a 75 kpc cometary 
emission of ionized gas trailing behind two Irregular galaxies. The Hq trails correspond in position and 
length with tails of syncrotron radiation. At the galaxy side opposite to the tails the two galaxies show 
bright HII regions aligned along arcs, where the star formation takes place at the prodigeous rate of 
~ IM yr~^. From the morphology of the galaxies and of the trailing material, we infer that the two 
galaxies are suffering from ram pressure due to their high velocity motion through the cluster IGM. We 
estimate that ~ 10^ M0 of gas, probably ionized in the giant HII regions, is swept out forming the tails. 
The tails cross each other at some 100 kpc from the present galaxy location, indicating that a major 
tidal event occurred some ~ 5 x 10^ yr ago. We exclude that mutual harassment produced the observed 
morphology and we show with numerical simulations that it could have marginally aided ram pressure 
stripping by loosening the potential well of the galaxies. 

Subject headings: galaxies: clusters: individual (A1367) — galaxies: irregular — galaxies: evolution — 
galaxies: intergalactic medium — methods: N-Body simulations 
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1. INTRODUCTION 

CGCG 97-073 and 97-079 (Zwicky et al, 1961-68) are 
Irregular members to the cluster of galaxies A1367. They 
have among the highest star formation rate found in Ir- 
regular galaxies: ^ IM yr~^ as derived from their H^ 
luminosity ~ 1.5 x lO^^erg sec~-^ (Kennicutt, Bothun & 
Schommer, 1984; Gavazzi et al. 1998). The star formation 
takes place in bright HII regions distributed along curved 
paths on the galaxy periphery facing the cluster center (see 
Fig. |l|). Both galaxies are significantly anemic on the op- 
posite side. Radio continuum observations with the VLA 
(Gavazzi & Jaffe 1985, Gavazzi & Jaffe 1987, Gavazzi et 
al. 1995), revealed the "head-tail" appearence of the radio 
source associated with these two galaxies. The tails ex- 
tend up to 75 kpc (assuming Hq = 75 km s^^ Mpc^^) on 
the side opposite to the bright HII regions. Observations 
in the 21 cm line of HI (Gavazzi 1989; Dickey & Gavazzi 
1991) revealed that both galaxies have a slightly deficient 
HI content, displaced in the direction marked by the ra- 
dio continuum tails, as opposed to their H2 content which 
appears normal in all respects (Boselli et al. 1994). These 
asymmetries suggest that the two galaxies are experienc- 
ing ram pressure due to their high velocity motion through 
the IGM. Abadi et al. (1999), Murakami & Babul (1999) 
and Quilis et al. (2000) performed high-resolution hydro- 
dynamical simulations of galaxies subject to ram pressure 
stripping in rich clusters. In less than 1 Gyr their galaxy 
models lose all the g result of ram pressure and 

viscous stripping when the density of the IGM and the 
transit velocity are as high as in the Coma cluster. Ex- 
tended gaseous tails form and the gas is shocked at the 

^Based on observations taken with the INT, WHT, NOT and SPM 
Spanish Observatorio del Roque de Los Muchachos of the I.A.C. The 



leading edge of the galaxies, considerably enhancing its 
density, thus most likely leading to intense star formation, 
as observed in 97-073 or in NGC 4522 (Kenney & Koop- 
mann 1999) in the Virgo cluster. 

These phenomena might contribute significantly to the en- 
richment of the intergalactic medium in clusters, which 
remains an unsettled issue (e.g. Madau, Ferrara & Rees 
2001; Mori, Ferrara & Madau, 2001; Silich, et al. 2001; 
Recchi et al. 2001). 

While taking deep (one hour exposure) Hq, images of 
the cluster A1367 with the Wide Field Camera at the 
Isaac Newton Telescope (La Palma), we serendipitously 
discovered a low surface brightness Hq, emission trailing 
behind 97-079, spatially coincident with the radio contin- 
uum tail. This emission was too faint to be measured 
and required confirmation. We thus took more observa- 
tions^ using narrow-band Hq filters: Ih with the NOT and 
4h with the SPM telescopes. Each individual observation 
confirmed the existence of the trailing Hq, . The resulting 
stacked 6h exposure image, which we present in this Let- 
ter, is sufficiently deep to allow a robust determination of 
the flux in the tail of 97-079 and revealed an even fainter 
trail behind 97-073. The cometary Hq trails discovered in 
A1367 are the most extended low-brightness Hq emission 
features ever detected. Other prominent examples are the 
the Magellanic stream (Weiner & Williams, 1996) and the 
cloud in Leo (Reynolds et al. 1986). 

2. He OBSERVATIONS 

We observed the field centered on CGCG 97-073 and 97- 
079 using three telescopes: the 2.5 m Isaac Newton Tele- 
telescopes. The INT, WHT and NOT are operated by the I.N.G. at the 
SPM telescope belongs to O.A.N de Mexico. 



2 



Table 1: The imaging instrumental set-up 



Telescope 




Date 




Pix 


Filter 


Tint 


seeing 










arcsec 


A 


sec 


arcsec 


INT 


26 


Apr 2000 


4 X 2048 X 4100 EEV 


0.33 


6725 (80) 


3600 


1.7 


INT 


26 


Apr 2000 


4 X 2048 X 4100 EEV 


0.33 


(Gunn) r' 


900 


1.7 


SPM 


23 


Apr 2001 


1024 X 1024 Thompson 


0.38 


6723 (80) 


14400 


1.9 


SPM 


23 


Apr 2001 


1024 X 1024 Thompson 


0.38 


(Johnson) R 


1800 


1.9 


NOT 


25 


Apr 2001 


2048 X 2048 Loral 


0.18 


6725 (60) 


3600 


1.1 


WHT 


9 


Fe6 2001 


2 X 2048 X 4100 EEV 


0.24 


6736 (48) 


6600 


0.85 


COMB 








0.38 


ON 


21600 


1.9 


COMB 








0.38 


OFF 


2700 


1.9 



scope (INT), the 2.5 m Northern Optical Telescope (NOT) 
at La Palma (Spain) and the 2.1 m telescope at San Pedro 
Martir (Mexico) (see Table 1). The observations were per- 
formed through narrow band filters centered at 6725 A, 
covering the redshifted Hq, and [Nil] lines. The underly- 
ing continuum was taken through broad band red filters 
(see Table 1 for details). The images were obtained in 
photometric conditions with a seeing ranging from 1.1 to 
1.9 arcsec. Each integration was split in shorter expo- 
sures (typically 20 min) to get rid of the cosmic rays. The 
photometric calibration was obtained exposing the spec- 
trophotometric star Feige 34. The individual images were 
bias subtracted and flat-fielded using combinations of ex- 
posures of several empty fields at twiglight. After back- 
ground subtraction, the images were combined. This was 
accomplished using IRAF tasks. Based on approximately 
100 common stars found in the field using DAOFIND, the 
proper coordinate transformations i.e. rotation, trasla- 
tion and re-sampling (independently on X and Y) were 
applyed using GEOMAP. The resulting stacked frames, 
corresponding to 6 hours and 45 min of integration time 
(ON and OFF band respectively), have 0.38" pixels. The 
intensity in the combined OFF-band frames was normal- 
ized to that of the combined ON-band one using the fiux 
ratio of several field stars. The combined NET-image was 
obtained by subtracting the normalized OFF-band frame 
from the ON-band one. The resulting 5' x 5' ON- and 
NET-frames are shown in Fig. 0. One additional Hq, ex- 
posure of the field was obtained with the 4.2m William 
Herschel Telescope (WHT) at La Palma. Given the excel- 
lent seeing conditions (0.85 arcsec) it provided high resolu- 
tion images of the two galaxies shown in Fig. ^ However, 
due to the bright moon conditions, it did not provide any 
useful information about the low-brightness features, thus 
it was not stacked with the remaining observations. 

3. Hc« TAIL PARAMETERS 

The Hq, + [Nil] fluxes within the optical extent of 
the two galaxies under study are found consistent with 
LogF(Ha-|- [Nil]) = -12.75 (-12.66) erg cm'^sec"! (here- 
after the first quantity refers to 97-079, the second one 
in parenthesis to 97-073) (Gavazzi et al. 1998). The 
flux estimate of the low brightness tails requires a care- 
full assessment of the quality of the flat-fielding. For ex- 
tended sources the dominant source of error is associated 



with the variations of the background on scales similar 
to the source. We measured the background in several 
10 X 10 arcsec^ regions around the field (comparable with 
the cross section of the tails) and determined the rms of 
the background on this scale. The flux density in the tails 
was found to decline from 8 tr to 2 cr at 3 ('^ 2) arcmin 
from their parent galaxies. Assuming a distance of 86 Mpc 
to A1367 the total extent of the tails results in 75 (50) 
kpc. The cross section of the tails is '--^ 8 kpc. The total 
flux in the trails, obtained integrating the counts above 
2 cr in rectangular regions of approximately 3 x 0.3 ar- 
cmin (2 X 0.3 arcmin) (excluding the galaxies themselves) 
results in: LogFtraii = —14.11 (-14.85) erg cm~^sec~^. 
The line intensity is 0.25 (0.07) Rayleigh (1 Rayleigh = 
10^/477 photons cm^^ s^^ sr^^), corresponding to an emis- 
sion measure (EM) of 0.67 (0.18) cm^® pc if the gas is op- 
tically thin at 10"* K. Assuming that the tails have cylin- 
drical symmetry, with a filling factor of 1, this implyes a 
mass of ionized gas of 9.6 x 10* (3.2 x 10*) M0, thus a 
plasma density of 9 x lO"'^ (5 x 10^'^) cm^^ or lower if 
the filling factor is < 1. This estimate is comparable with 
the present content of neutral hydrogen of the two galax- 
ies: 1.6 X 10^ (2 X 10^) M0 respectively (Gavazzi, 1989). 
The plasma in the tails was probably ionized inside the 
galaxies by the actively star forming HII regions on their 
leading-edges. However it must have originally consisted 
of neutral hydrogen, therefore a rough estimate of the total 
gas loss from the two galaxies can be derived using the HI 
deficiency parameter, as defined by Giovanelli & Haynes 
(1985). Although one should take this estimate cautiously, 
because for Irregular galaxies the reference HI content is 
highly uncertain, from the HI deficiency parameter 0.25 
(0.16) of the two galaxies respectively we estimate that 
1.25 X 10^ (9.1 X 10*) M0 of gas was left behind in the 
tails or that the two galaxies have lost approximately 40 
% of their original gas. This estimate is consistent with 
the mass previously found in the ionized gas. 
Assuming a transit velocity through the cluster of ~ 
1200 km s~^, as derived from V2 x avei (here we assume 
that the cluster can be modeled as an isothermal sphere, 
as in Cayatte et al. 1994) with a-yei — 822 km s^^ for 
A1367 (Struble & Rood 1991), from the length of the 
tails we derive that the ionized material survived some 
10'' *(10^'®) yr. The recombination time = l/N^aA, 
where aA — 4.2 x 10"-^^ cm^s~^ (Osterbrock, 1989) is 
10^ °(10^'^) yr i.e. about 5 times shorter than the survival 
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time. The two times would however become eonsistent if 
the plasma density was about 5 times lower than the value 
estimated above, depending on the filling factor. Although 
it is likely that the gas in the tails came out ionized from 
the galaxies, it cannot be excluded that the eddies along 
the tails contain sufficient turbulent energy to currently 
sustain its ionization along the tail. 

4. DISCUSSION 

Both 97-073 and 97-079 appear to have lost to the tails 
^ 40% of their original gas content. Stripping of the gas 
component of galaxies in clusters can occur either due 
to the ram pressure exerted by the intracluster medium 
(Gunn & Gott 1972) or because of tidal interactions. Ram 
pressure acts on a very short timescale, of the order of lO'' 
yr, and produces a stream of gas that trails the galaxy 
(Abadi et al. 1999; Quilis et al. 2000). Tidal interac- 
tions within galaxy clusters are of two different types; (a) 
galaxies interact with the global potential of the cluster, 
being tidally shocked each time they approach the peri- 
center of their orbit inside the cluster core (Merritt, 1983), 
and (b) they interact with the other galaxies through re- 
peated high-speed fly-by encounters, a mechanism known 
as galaxy harassment (Moore et al. 1996, 1998, 1999). 
Global tidal interactions can remove a large fraction of the 
gas reservoir of galaxies (Mayer et al. 2001a, b) but they 
require a long timescale, of the order of 7-10 Gyr, which 
is longer than the typical age of galaxy clusters (Rosati 
et al. 2000). Galaxy-harassment acts on a much shorter 
timescale in clusters; the collision rate per galaxy in a rich 
cluster like Coma is ^ 1/Gyr and 1-2 collisions with a 
large galaxy might remove ^ 50% of the gas in a disk 
galaxy (Moore et al. 1998). However, harassment, as any 
other tidal mechanism, is expected to produce a leading 
gaseous tail in addition to the trailing tail, which can be 
excluded for the galaxies here analyzed. Overall, ram pres- 
sure seems to provide the most likely explanation for the 
features observed. 

Yet, given the density of the gas in the NW subclus- 
ter of A1367, namely ~ 7 x 10""* atoms cm~^ (Don- 
nelly et al. 1998), a factor ~ 5 lower than the density in 
the core of Coma, ram pressure is expected to be rather 
weak according to recent numerical simulations (Abadi et 
al. 1999; Quilis et al. 2000). However, these studies 
focused on the effects on large spiral galaxies whose po- 
tentials are much more centrally concentrated and should 
yield stronger restoring forces than those of the faint, ir- 
regular galaxies considered here. Moreover, the relative 
positions and present apparent direction of motion of the 
galaxies (we assume that the galaxies are moving in the 



opposite direction with respect to the gaseous tails) are 
consistent with them having undergone a very close en- 
counter; assuming that their velocity is ~ 1200 km s~^, 
the encounter likely took place some 10^'^ yr ago. The 
tidal encounter could have lowered the restoring force by 
loosening the potential well of the galaxies; we tested this 
hypothesis with a few N-Body simulations performed with 
the binary treecode PKDGRAV (Dikaiakos & Stadel 1996; 
Stadel & Quinn, in preparation). The target galaxy is rep- 
resented by a high-resolution N-body model analogous to 
those described in Mayer et al. (2001a, b), that obey the 
TuUy-Fisher relation (Zwaan et al. 1995), and comprises 
an exponential stellar disk of 50,000 particles embedded in 
an isothermal dark matter halo of 400,000 particles. The 
model has structural parameters and a luminosity con- 
sistent with those measured in 97-073 and 97-079. The 
perturbing galaxy is represented by the external poten- 
tial of an isothermal halo with mass comparable to that 
of the target model. The relative velocity of the galax- 
ies is varied in the range Vrei = 500 — 1000 km s^^ (i.e. 
comparable to the velocity dispersion measured for the 
galaxies in the cluster, see section 3), the impact parame- 
ter is equal to the core radius of the perturber and various 
disk orientations are considered. We analyze the target 
galaxy ^ 10^ yr after the collision in the various runs and 
find that its mass distribution, and thus its restoring force, 
has been only slightly affected, both because the pertur- 
bations are intrinsically small and because the response of 
the system is slow compared to the time elapsed since the 
collision has occured. However, if the galaxies travel at 
~ 1200 km/s in the cluster, the ram pressure would over- 
come their restoring force at more than the two disk-scale 
lengths from the center, mainly because their potential 
is intrisically shallow. Stripping down to 2 scale lengths 
would remove ^ 25% of the gas if it were distributed as 
the stellar disk; turbulence and viscosity would enhance 
gas stripping (Quilis et al. 2000) and, in addition, gaseous 
disks in late-type galaxies arc usually more extended than 
the stellar disk of our model (de Blok & McGaugh 1997) 
and would be more fragile. Hence, unless 97-073 and 97- 
079 have potentials much more concentrated than assumed 
here, our results are consistent with the hypothesis that 
ram pressure has produced their gaseous tails. 

G.G. wishes to thank Christian Bonfanti for his contri- 
bution to the image reduction and Luca Cortese for use- 
ful discussions. L.M. thanks Thomas Quinn and Joachim 
Stadel for allowing him to use PKDGRAV and Ben Moore 
for useful discussions. L. Carrasco research is supported 
by CONACYT research grant G28586-E. 
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Fig. 1.— The + [Nil] NET image taken with the WHT with 0.85 arcsec seeing showing the bright parts of 97-079 
(left) and 97-073 (right). (J2000) celestial coordinates are given. The center of the cluster is at S-E. An arrow marks the 
direction of the low brightness tail (see Fig. 2). 
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Fig. 2. — The stacked 6 hours Hq, + [Nil] ON band (left) and NET exposures (right) of the region containing the two 
galaxies under study shown at high contrast to enhance the extended tails. The upper corners of the images suffer from 
filter vignetting. 



